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Abstract 

We consider the problem of performing the preliminary "symmetry 
classification" of a class of quasi-linear PDE's containing one or more 
arbitrary functions: we provide an easy condition involving these func- 
tions in order that nontrivial Lie point symmetries be admitted, and 
a "geometrical" characterization of the relevant system of equations 
determining these symmetries. Two detailed examples will elucidate 
the idea and the procedure: the first one concerns a nonlinear Laplace- 
type equation, the second a generalization of an equation (the Grad- 
Schliiter-Shafranov equation) which is used in magnetohydrodynamics. 



Key words: symmetry classification; quasi-linear PDE's; symmetry deter- 
mining equations; generalized Laplace equation; Grad-Schluter-Shafranov 
equation 



Running title: Symmetry Classification of quasi-linear PDE's 



1 



1 Introduction 



The analysis of symmetry properties of differential equations is a well estab- 
lished and widely used tool both for studying general properties of the equa- 
tions and for finding their solutions (see e.g. PQ-[I] and references therein); 
actually, the determination of Lie-point symmetries (we will consider only 
symmetries of this type) is by now an almost completely standard routine, 
thanks also to some suitably dedicated computer packages (see e.g. [8]- [10j) . 

The situation is however considerably different when one deals with the 
problem of performing the "symmetry classification" of an equation which 
contains one or more arbitrary functions, and one wants to discover how the 
symmetry properties depend on the choice of these functions: the discussion 
may be far from easy, as far as the symmetries may drastically change when 
the functions are changed: see e.g. |llj-|15j and references therein. 

In this paper, we discuss the case of PDE's of the form given below 
(eq. ([!])) containing one or more arbitrary functions F^{u) of the unknown 
variable u, and we shall give an easy condition involving these functions 
in order that the equation admits nontrivial symmetries. This will provide 
a neat "geometrical" characterization of the relevant system of equations 
determining these symmetries, and a direct way to perform the complete 
symmetry classification of the given equation. Two examples will elucidate 
the idea and the procedure. The first one concerns a nonlinear Laplace-type 
equation; the second example deals with a generalization of an equation (the 
Grad-Schliiter-Shafranov equation) which is used in magnetohydrodynamics 
and plasma physics |16j . 

2 Preliminary results 

For the sake of simplicity we consider only second order equations for the 
unknown function u = u(x, y) of two independent variables x, y (but the 
extension to more general cases is completely straightforward), and we will 
deal with quasi-linear PDE's of the following form 
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Oil Uxx + 0,12 U xy + a 2 2 Uyy + h u x + b 2 u. 



■v — 





or, in a short-hand notation, 



£[u] = aeF e (u) 
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where = aij(x,y), hi = bi(x,y), ct£ = otg(x,y) are given (smooth) func- 
tions, and Fi{u) are L arbitrary (smooth) functions of u (in the examples 
below we will deal with just one or two functions Fg(u)). It is understood 
that no linear relations exist between ag and between Fg. We will exclude 
from our consideration the rather trivial case when Fg{u) are linear func- 
tions of u, which usually can be more simply and conveniently discussed 
separately by means of direct calculations. In this paper we are looking 
only for Lie-point symmetries of ([I]); we will denote their Lie generator by 

X = £( x ,y,u)-^ + 7]{x,y,u)^- + ip(x,y,u)-^ . 

Assume e.g. an ^ (and put then an = 1). 

Although essentially standard (see [TJ [2j [7j , and also [17] for some gen- 
eralization), let us summarize for convenience and in order to fix notations, 
these first basic results, which can be easily obtained imposing the usual 
symmetry condition 

X( 2 )(A) = 
A=0 

where A = £[u] — a^Fg, and X^ is the second prolongation of X. 

Lemma 1 For any choice of the functions Fg, the coefficients of the Lie- 
point symmetry operators admitted by a PDE of the form ([I]) satisfy the 
conditions £ u = t] u = 0; (p uu = 0, i.e. 

Z = €(x,y) , rj = r](x,y) , ip = A{x,y) +uB(x,y) (2) 

and the unique determining equation which involves the functions Ft{u) takes 
the form (with F[ = dFg/du, sum over £ = 1, . . . , L) 

Pe(x, y)F t (u) + p L +e(x, y)F' t (u) + p 2 L+e(x, y)uF' l (u)+ (3) 
P3L+i{x, y)u + P3L+2(x, y) = 

where the coefficient functions Pi(x, y) (i = 1, . . . , 3L + 2) are given by 

Pi = Ba e - (£ 2 + r] y )a e - aniyOt-t - ia^ x - r) a e>y , (4) 

PL+l = -at A , P2L+1 = -OL(B , P3L+1 =£[B) , P3L+2 = £[A] 

with £ x = d^/dx, OLi^ x = dai/dx etc. (£ = !,... ,L). 
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Considering now the determining equation ([3]), and observing that the pi 
depend only on x, y and not on u, one immediately realizes that, if the 
3L + 2 functions /j defined by 

f = (F e ,F^,uF^u,l) (£ = 1,...,L) (5) 

are linearly independent, then ([3]) can be satisfied if and only if 

Pi = (i = l,...,3L + 2) . (6) 

Recalling now the definition of kernel of the full (or principal) symmetry 
groups pQ of eq. JT]), i.e. the intersection of all symmetry groups admitted 
by ([1]) for any arbitrary choice of Fg(u), we can then state the following 
property. 

Lemma 2 Conditions (El) characterize the kernel of the symmetry groups 
of equation |7J). 

Indeed, conditions ©, together with the other determining equations (not 
involving Fi), determine the functions ^,rj,A,B (i.e. the symmetries ad- 
mitted by equation ([1])), which are independent of the the choice of the 
functions F^. These symmetries may be considered "trivial" in this context: 
for instance, if all coefficients a, b, a in ([1]) are independent of y, then such 
a symmetry operator is d/dy. Some not so obvious examples of symmetries 
of this type will be presented later. 

Therefore, a first conclusion is that, in order to have "nontrivial" sym- 
metries (i.e. really dependent on the choice of the functions Fg ), a necessary 
condition is the existence of some linear dependence among the functions ([5]) . 

Another relevant remark which will emerge from our discussion is the 
important role played also by the coefficient functions ai(x,y) in the deter- 
mination of the admitted symmetries. 

3 Conditions for the existence of symmetries. 

Consider the linear space generated by the 3L + 2 functions fi defined in ([5]) , 
and, according to our above remarks, now assume that there are some linear 
relations among these functions. Then the fi span a space with dimension 
k < D = 3L + 2; if this is the case, the D coefficients pi are forced, accord- 
ing to ([3]), to belong to the orthogonal (D — A:)-dimensional subspace (with 
respect to the standard scalar product in R/°), and the functions Pi(x,y) 



4 



turn out to be subjected to k linear conditions. For instance, if there is just 
one linear relationship between the /«, say 

D 

= (7) 

i=i 

where not all the constants Aj vanish, then k = D — 1 and the functions pi 
span a 1-dimensional subspace and must satisfy D — l equations of the form 
(assuming that, e.g., Ad ^ 0) 

pi\ D = \ip D , p 2 Xd = Mpd , ■■■ , PD-l^D = ^D-lPD • (8) 

We can then state our main conclusion, which characterizes the crucial de- 
termining equation which contains the functions in the following "geo- 
metrical" form. 

Proposition 1 Eq. (QP admits nontrivial symmetries only if the D = 3L+2 
functions {5]J are linearly dependent. If this is the case, the D functions 
Pi{x,y) given by appearing in the determining equation f3|) span the 
subspace orthogonal (with respect to the standard scalar product in H D ) to 
the k— dimensional (k < D) subspace spanned by the functions The 
admitted symmetries are completely determined by imposing this orthogo- 
nality condition to the coefficient functions Pi(x,y), together with the other 
determining equations not involving the functions Fg(u). 

Before considering explicit examples, let us remark that the complete 
symmetry classification must be accompanied by the determination of the 
equivalence group i.e. the group of the transformations which leave 
invariant the differential structure of the PDE. Standard calculations show 
easily that, for any fixed choice of the functions a^-, bi, ct£ in equation ((!]), 
the equivalence group includes in particular, expectedly, the scalings u — > 
cm, Fi — ► cFi and the translation u — ► u + Co (c, Co = const.). Other 
transformations involving also the variables x, y can appear for particular 
choices of the functions a^, bi, ot£. The transformations belonging to the 
equivalence group will play an important role in performing the complete 
symmetry classification of our equations. 

4 First example: a generalized Laplace equation 

To illustrate the main idea and the procedure, and also to clarify some 
details, we are now going to examine some examples, which can be notewor- 
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thy also for their different and interesting peculiarities. We start considering 
the simplest case, where the r.h.s. of ([T]) contains only one term a(x, y)F(u) 
(then D = 5). 

First of all, let us remark that, independent of the form of the equation 
(PQ), not all linear combinations among the /j defined in ([5|) can be chosen 
arbitrarily. For instance, no relation exists between u and 1, and also, 
having excluded the trivial case of linear F(u), between F, u, 1. On the 
other hand, the necessary linear dependence between the functions fi implies 
immediately that the presence of some symmetry is possible only if F(u) is 
an exponential or a power of u. It can also happen that more than one linear 
relation holds: e.g., if F = u 2 /2 + u, then both relations 

F'-u-l = and 2F-uF'-u = (9) 

hold simultaneously, i.e. k = 3. In this case, the orthogonality condition 
expressed by Proposition 1 becomes 

Pi + 2p 4 - 2p 5 = , p 2 + P5 = , p 3 - Pa + P5 = . (10) 

As an explicit example, let us consider the following generalization of 
the classical nonlinear Laplace equation 

A = V 2 -u - a(x, y)F(u) = (11) 

where a(x, y) is a given function. In this case, the determining equations 
not containing F imply in particular 

ix = % , i y = -r]x , B = const (12) 

whereas the crucial determining equation ([3j) involving F is 

p x F + p 2 F' + P3 uF' +P4U + P5 = (13) 

with the coefficient functions pi(x,y) given by 

Pi = Ba- (Cx + %)« - (£,a x + r]a y ) , p 2 = -aA , (14) 
p 3 = - a B , pi = V 2 B = , p 5 = V 2 A . 

Let us first discuss the kernel group. The conditions Pi = (i = 1, . . . , 5) 
characterizing the transformations in the kernel group (see Lemma 2) give 
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now A = B = tp = and the condition p\ = 0, which now reads ot(£ x + r] y ) + 
(£, a x + V a y) = 0- Introducing a harmonic function <I> = &(x,y) such that 

$x = Z ®y = -V 

this condition can be more conveniently transformed into an equation for 
the single unknown <I>: 

a(<$> xx - $yy) + a x $ x - a y <S> y = . (15) 

Solutions of this equation clearly depend on the choice of the function 
a(x,y). For instance, if a = const., it gives £ x + r\ y = 0, which, together 
with (|12p. implies that the symmetries in the kernel group are, as expected, 
only translations and rotations of the variables x,y. If a = exp(2x) then 
<3? = exp(— x)(pi cosy + C2siny) + c%y + C4 and the kernel contains, apart 
from the translation generated by d/dy, the transformations generated by 

/ d d \ ( d d 

Xi = exp(-x)(^cosy — -siny— J , X 2 = exp(-x)(^siny — +cosy— 

With a = x r , one has that if r 7^ — 2 then the kernel contains only the trans- 
lation generator d/dy, whereas if r = —2 it also contains the transformations 
generated by the two operators 

Xi = 2xy- (x 2 - y 2 )^- and X 2 = x— + y— . 

ox oy ox oy 

The first transformation describes the kernel group even if a = x~ 2 (5 (y / (x 2 + 
y 2 )), where (5 is an arbitrary function. 

It can be remarked, incidentally, that if we reverse the argument for a 
moment, one has that: given any harmonic function $ (and then any couple 
of harmonic conjugate functions £, 77), there are some a(x,y) which solve 
equation (|15p . and then, with these functions a, the kernel group contains 
precisely the symmetry generated by the operator X = ^(d/dx) + n(d/dy). 

Let us now finally perform the symmetry classification of eq. (jlip . As 
already remarked, its equivalence group may contain, in addition to the 
transformations listed at the end of Section 3, and depending on the spe- 
cific choice of the function a, other transformations possibly involving also 
x, y. As we shall see, however, these are not relevant for the symmetry 
classification of equation (jlip . 

According to our procedure, it is immediately seen that just one linear 
relation between the five fi can exist. For instance, in the case F = u 2 /2 + u 
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mentioned above, admitting the two linear relations Q, conditions (jlOp 
would lead to pi = 0, i.e. only the kernel symmetries. We then assume the 
existence of a single linear relation: 

XiF + X 2 F' + X 3 uF' + X 4 u + X 5 = (16) 

with not all Xi equal to zero. Observing that p 4 = 0, one has from ([8]) that 
A4 = 0. We now distinguish the cases A3 7^ and A3 = 0. 

Let A3 7^ 0. It is not restrictive to put A3 = 1, and (up to a translation 
of u) X2 = 0, which implies X5P2 = 0. If A5 7^ 0, then p 2 = would imply 
A = and also P5 = = P3A5; now, if ps = it remains only p\ 7^ 0, 
and from A5P1 = Aips = one concludes that A5 = 0. Therefore, we get 
X 1 F + uF' = 0, i.e. 

F(u) = u m with m = -Ai (m 0, 1) . 

Using now flSJ), which become p\ + mp^ = , p 2 = Pa = P5 = 0, we get 

A = and aB(m — 1) + + 77^)0 + (£0^ + r/a^) = . 

The last equation relates the symmetry coefficients £, rj, B with the specific 
form of the function a(x, y). If for instance a = x r , then £>(m — 1) + £ x + 
r/y = 0, but i? must be ^ 0, otherwise also p% = pi = 0, i.e. the kernel 
group. Therefore, £ must be proportional to x and equation (jXTJ) admits the 
symmetry operator 

x = ^- l) ( x Tx +y ^)- {r + 2)u lL 

(and obviously the translation of the variable y, and also the translation of 
x and the rotations of x, y in the case r = 0, i.e. if a = const.). 

Let now A3 = 0. Then necessarily A2 7^ 0, and one can put A2 = 1 and 
also Ai = 1 (possibly up to a scaling of u). Assume first A5 = 0; therefore, 
from (fT6|) . 

F(u) = exp(— u) 
and the conditions (JSj) , (fH|) for the functions pi become now 

P3 = —aB = , P5 = V 2 A = , and p\ = P2 i.e. 

{ix + Vy)a + {£a x + 7ja y ) = aA . 
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As before, we can consider some examples. If a = const., the last equation 
shows A = £ x + r\ y and then the most general symmetry of the equation 
V 2 u = exp(— u) is 

x = ^ + ^ + & + n«)l 

where £, rj are arbitrary harmonic conjugate functions: this is the well 
known case of the classical Liouville equation (in its "elliptic" form; for 
a full discussion of the symmetry properties and other related features of 
the Liouville-type equations, see [18j-[23j). If instead for instance a = y r 
then the admitted symmetry operators are 

X 1 = x^- + y^- + (2 + r)^- , X 2 = (x 2 -y 2 )^- + 2xy^- + (2r + 4)x-^- 
ox oy ou ox Oy Ou 

and the translation d/dx. Assume now A5 ^ 0, then 

F(u) = exp(— u) — A5 . 

Introducing the transformation u — > u + u, where u = u(x, y) satisfies the 
equation V 2 n + Asa = 0, one obtains the new equation 

V 2 u — 5exp(— u) = 

where a(x, y) := a exp(— u), which has precisely the same form as the equa- 
tion considered before. Without repeating details, it can be interesting to 
provide just one illustrative example. Let 

V 2 n - /cx" 2 (exp(-n) + l) = 

where k = const. It is easy to see that if k = 2 this equation admits the 
symmetries 

, d d r. 2 \ d 

x = £7r + r ?7r + & + %--£br 

ox oy V x / Ou 

where ^, 77 are arbitrary harmonic conjugate functions; if instead k 7^ 2, the 
admitted symmetries are only those in the kernel group. 

The above results concerning eq. (jlip can be stated in a complete form 
as follows. 

Proposition 2 Given a function a = a(x,y), consider this equation for the 
harmonic function $ = y) 

a($ xx - $ yy ) + a x $ x - a y $ y - aC = . (17) 
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Let £ = $ x , r] = —$> y . Assume first C = 0: for any solution of 
the kernel group of the generalized Laplace equation ill]) is generated by the 
symmetry operator 

dx ^ dy 

If F(u) = u m , for any solution of ( f7?| j with C = const.^ 0, eq. 177]) 
admits the symmetry operator 



If F{u) = exp(— u), for any solution 3> of (H\ ) with C = C(x,y) any har- 
monic function, eq. ( ti 1\) admits the symmetry operator 

d d d 

X = ^TT + ^TT + ^^'^TT • 
ox ay ou 

In particular, if a =const. then C = £ x + rj y where £, r/ are arbitrary har- 
monic conjugate functions, and the case of the standard "elliptic" Liouville 
equation is recovered. If finally F{u) = exp(— u) + c (c = const.), the above 
case is recovered by means of the transformation u — > u+u, where u = u(x, y) 
satisfies the equation V 2 u — ca = 0. This completes the symmetry classifi- 
cation of the PDE {HP , apart from the transformations in the equivalence 
group. 



5 An example with two arbitrary functions. 

We now consider the case of a PDE of the form (JTJ) with two arbitrary 
functions F^{u), i.e. L = 2, D = 8. To avoid excessive generality, let us 
restrict our study to a PDE of the following form 

u xx + Uyy + \ x = a(x,y)Fi(u) + F 2 (u) (18) 

here b\ = a/x, a ^ is a constant, a\ = a(x, y) a given function and «2 = 1- 
The choice of this equation is motivated and suggested by the theory of 
plasma physics: it is indeed a generalization of the Grad-Schluter-Shafranov 
equation (see [K]), which is obtained putting in (fT8l) a = — 1, a = x 2 , 
and describes the magnetohydrodynamic force balance in a magnetically 
confined toroidal plasma. In this context, u is the so-called magnetic flux 
variable, x is a radial variable (then x > 0), while the two arbitrary functions 
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-Fi(u), F2{u) are flux functions related to the plasma pressure and current 
density profiles. 

The determining equations not involving Fg(u) give in this case 

ix = Vy, iy = -Vx, and B = + b, b = const . 

First of all, the kernel group is immediately seen to be trivial (apart obvi- 
ously from the translation generated by d/dy, in the case where a depends 
only on x; the possible presence of this symmetry will be tacitly understood 
in the following). Indeed, from pi = (see Lemma 2), one has A = B = 0, 
then the above equation implies £ = (2b/ a) x, and the condition p 2 = with 
a 2 = 1 gives finally £ = 0. 

Let us now start assuming that there are exactly two linear relations 
involving the functions F\ and F<i separately: 

AiFi + A 3 F{ + X 5 uF{ + X 7 u + A 8 = (19) 
A 2 F 2 + A 4 i^ + XquF^ + fi 7 u + H8 = 

with not all \i,fJ,i equal to zero. 

Let A5A6 7^ 0, and put A5 = Xq = 1. According to Proposition 1, the 
symmetry coefficients pi, given by (J2J), satisfy then the six linear conditions 

Pi = A1P5 , P2 = A 2 £>6 , P3 = A3P5 , Pa = A4P6 , (20) 

Pi = A7P5 + MPG , P8 = A 8 P5 + 

With A5 7^ 0, we can put A3 = 0, up to a translation of u. Conditions ([20]) 
and the expression of the coefficients pi give ps = 0, A = 0, then p^ = pg = 
and therefore also A4 = Xg = /is = (indeed, p^pe ^ 0, otherwise all pi = 0); 
condition p 2 = A 2 j>6 implies that £(x,y) must satisfy an equation of the form 

6r = &o — \- fei (koi ki = const) 

x 

which admits harmonic solution only of the form £ = cx, c = const. On 
the other hand, Ai = P1/.P5, A 2 = P2/P6 imply that a is forced to satisfy 

xa x + ycty 

— ? = r = const . 21 

a 

This means that if a does not satisfy this condition, no symmetry is allowed; 
we then assume for a the form 

a(x,y) = x r (3{y/x) 
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where j3 is arbitrary. Notice that, with a of this form, a new transformation 
is included in the equivalence group, namely the scaling x — > cx, y — > 
cy, F\ — > c 2 ~ r Fi, F 2 — > c~ 2 F 2 . We also deduce i? = const. ^ 0, p 7 = 0, 
which implies in turn Xj = fij = 0. Then we are left with 

Ai.Fi + uF[ = , X 2 F 2 + uF' 2 = 

giving (thanks to some scalings - all these transformations belong indeed to 
the equivalence group) F\ = u~ Xl , F 2 = u~^ 2 where 

-Ai = 1 - -^(2 + r) = 1 + , -A 2 = l + - with B = —cq 

B q q 

with admitted symmetry generated by 

d d d 

X = x a~ + ya~ ~Q U TT ■ 
ox ay ou 

Let now A5 = Xq = 0. Then necessarily A3A4 7^ 0. According to Propo- 
sition 1, the orthogonality condition now reads (with A3 = A4 = 1) 

Pi = Aip 3 , p 2 = X 2 p4 , p 5 = p G = , p 7 = A 7 p 3 +^7P4 , p 8 = A 8 p 3 +/i 8 p4 • 

In this case, one has immediately B = 0, then pj = and A7 = ^7 = 0, 
and again ^ = cx. From p 2 = p^ one has 2£ x = A = const., which gives 
Ps = A§ = fig = 0. Up to a scaling of u, one can choose X 2 = 1, the equations 
for F\, F 2 are then 

A1F1 + F[ = , F 2 + F' 2 = 

giving F% = exp ( — Xiu) , F 2 = exp ( — u), and finally from Ai = p\/pz 
one deduces the same condition ([2"T]) as before for the function a(x,y), and 
Ai = 1 + (r/2). 

The conclusion will be stated in complete form in the following Proposi- 
tion 3. Indeed, with the same, and even simpler, arguments used in the two 
above cases, it is an easy task to see that no other possibilities are left to 
the PDE (|18p of admitting symmetries. 

Proposition 3 The kernel group of equation U8\) is trivial ( apart from the 
translation y — > y + c if a depends only on x). Except for this, a necessary 
condition in order that equation 17%)) may admit symmetries is that the func- 
tion a has the form a(x,y) = x r (3{y/x), where (3 is an arbitrary function. 
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With a of this form, equation \18\) admits a symmetry only with the fol- 
lowing choices for the functions F\ (u) , F 2 (u) ( up to transition to equivalent 
functions via equivalence group): 

a) Fi(«) = u 1+ ( r+2 )/« , F 2 (u) = u 1+ ^ 

for all q ^ 0, with admitted symmetry operator 

d d d 

x = + y* — i u ^~ 

ox ay ou 



and 

b) F 1 (u) = exp (- (1 + ^) uj , F 2 {u) = exp(-n) 

with symmetry operator 



d d d 
X = x— + y—+2— . 

ox oy ou 

As stated at the beginning of Section 2, we have excluded from our analysis 
the case where the functions F± and F% are linear functions of u. This case 
indeed, in general of minor interest, can be more conveniently considered 
by means of a separate and direct calculation. Specifically, if the r.h.s. of 
equation ([18]) has the form a(x,y)u + 0(x,y), the admitted symmetry is, 
not surprisingly, generated by 

X = (cx + f(x,t/))^ 

where c is a constant and *$>(x,y) is any solution of the PDE 

= 5* - C f3 . 

In the particular case of the Grad-Schluter-Shafranov equation of plasma 
physics |16j . the above results were already presented in Ref. [23], but with- 
out any details in the calculations and without any reference to the proce- 
dure, which is instead one of the main purposes of the present paper. In 
the same reference [23] one can also find some physical comment on the 
symmetry properties of the above equation. 
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